Low energy states dynamic of entanglement for spin systems * 
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We have composed the ideas of quantum renormalization group and quantum information by 
exploring the low energy states dynamic of entanglement resources of a system close to its quantum 
critical point. We demonstrate the low energy states dynamical quantities of the one dimensional 
magnetic systems could show the quantum phase transition point and shows the scaling behavior 
in the vicinity of the transition point. To present our idea, we study the evolution of two spins 
entanglement in the one-dimensional Ising model in the transverse field. The system is initialized 
as the so-called thermal ground state of the pure Ising model. We investigate evolvement of the 
generation of entanglement with increasing the magnetic field. We have obtained that the derivative 
of the time at which the entanglement reaches its maximums with respect to the transverse field, 
diverges at the critical point and its scaling behaviors versus the size of the system are as same as 
the static ground state entanglement of the the system. 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 

A fundamental difference between quantum and clas- 
sical physics is the possible existence of nonclassical cor- 
relations in quantum systems. The physical property re- 
sponsible for this quantum correlation is called Entan- 
glement Entanglement has been recognized as an 
important resource for quantum information and com- 
putation However the role of entanglement in quan- 
tum phase transition (QPT) [3| is of considerable inter- 
est. QPT as well as classical ones are characterized by 
detecting nonanalytic behaviors in some physical proper- 
ties of the system. It is often accompanied by divergence 
in some correlation functions, but quantum systems pos- 
sess additional correlations which do not exist in a clas- 
sical counterpart, the entanglement. Entanglement is a 
direct measure of quantum correlations and shows non- 
analytic behavior such as discontinuity in the vicinity 
of the quantum phase transition point H[. An im- 
portant motivation to study the interconnection between 
condensed matter and quantum information is to inves- 
tigate whether it is possible to better characterize con- 
densed matter states by looking at their entanglement 
properties. Recently, there has been extensive analysis 
of entanglement in quantum spin models |6j. Various 
models were considered for entanglement generation and 
their static Q as well as dynamical properties Qwere 
investigated. A thorough understanding of the dynam- 
ical evolution of entanglement in the spin models has 
obviously implications for the performance of quantum 
information processing, as well as for understanding of 
fundamental quantum mechanics. 

Our main purpose in this work is to compose the 



ideas of quantum renormalization group Q and quan- 
tum information theory to study the evolution of the 
dynamical properties of the spin models in low energy 
states. To have a concrete discussion, the one dimen- 
sional S = \ Ising model in transverse field (ITF) has 
been considered by implementingthe quantum renormal- 
ization group (QRG) approach |9l-[llj. 



II. QUANTUM RENORMALIZATION GROUP 

The main idea of the RG method is the mode elimina- 
tion or thinning of the degrees of freedom followed by an 
iteration which reduces the number of variables step by 
step until reaching a fixed point. In KadanofFs approach, 
the lattice is divided into blocks. Each block is treated 
independently to build the projection operator onto the 
lower energy subspace. The projection of the inter-block 
interaction is mapped to an effective Hamiltonian (H e ^) 
which acts on the renormalized subspace [T3 - 4T3 | . 

We have considered the ITF model on a periodic chain 
of N sites with Hamiltonian 
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where J > is the exchange coupling and g is the 
transverse field. From the exact solution [la ] it is known 
that a second order phase transition occurs for g c = 1 
where the behavior of the order parameter or magneti- 
zation is given by < a x >= (1 — g) 1 ^ 2 for g < 1 and 
< g x >= for g > 1. 

To implement QRG the Hamiltonian is divided 
to two-site blocks, H B - ^ N/2 uB — ^ uB 
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n the inter-block part, H BB = 
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a-component of the Pauli matrix at site j of the block 
labeled by /. The Hamiltonian of each block (hf) is 
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diagonalized exactly and the projection operator (Pq) 
is constructed from the two lowest eigenstates, Pq = 
I'^oXV'ol + where |-0o) is the ground state and 

\tjji) is the first excited state. In this respect the effec- 
tive Hamiltonian (H e ^ = P [H B + H BB ]P ) is similar 
to the original one (Eq.([T])) replacing the couplings with 
the following renormalized coupling constants. 



J' = J- 
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III. REDUCED DENSITY MATRIX AND 
EVOLUTION OF CONCURRENCE 

The entanglement is a local quantity which includes 
the global properties of a system. Generally, the global 
properties of a system enters the entanglement effectively 
by summing over the whole degrees of freedom except the 
local one. In other words, a system can be supposed of 
a single site and a heat bath (the rest of system). It is 
supposed that the effect of a heat bath can be replaced by 
an effective single site quantity, the entanglement. The 
effective single site represents the long range properties of 
the model and not the microscopic ones. Therefor we can 
enter the global properties of the model to entanglement 
(the local quantity) using the renormalization group idea. 
In this respect, we always think of a two site model which 
can be treated exactly. However, the coupling constants 
of the two site model are the effective ones which are 
given by the renormalization group procedure. This can 
be used as an new method to calculate the low energy 
states dynamic of entanglement in a large system. 

The two site Hamiltonian of ITF model in the space 
spanned by {| ft), I U), Ut>, ] U)} (I t) and | |) denote 
the eigenstates of er 2 ), can be expressed as 
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FIG. 1: (Color online) Evolution of the concurrence under 
RG versus g for t = 1.5 and in the inset for t — 0.4. 



The density matrix for the two sites system at time 
t is p(t) — U{t)p(Q)U\ where p(0) is the the density 
matrix of system at t = which is the thermal equi- 
librium of the system. We choose the thermal ground 
sate of the pure Ising model (no transverse field applied) 
p(0) = ±(\a x ,+)(cr x ,+\ + \a x ,-)(a x ,-\) as the initial 
state which the concurrence equals zero. 

Here \a x , +) and \a x , — ) are the two degenerate ground 
state of Ising model with all spins pointing respectively 
to the positive and negative x direction. 

For bipartite entanglement, a commonly used measure 
for arbitrary state of two qubits is the so called concur- 
rence 03. The concurrence is defined as 
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However the time evolution operator U(t) — e of 
the two site Hamiltonian (Eq.Q) has the following form 



U(t) = J 
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where 



, 29 sin( Jt y/l + 4g 2 ) + cos( Ji \/l + V) 

V 1 + V 



: sin 



u m (t) = u 41 (t) = -= t 
vi + V 

U 22 {t) = U 33 (t) = cos J t, U 23 (t) = U 32 (t) =isinJt, 

U AA {t) = u* v 



C{t) = max{0, 2X max (t) - tr^J p{t)p{t)} 
where p(t) = (a y ®a v )p*{t){a v (gfcr 2 ') 



(5) 



largest eigenvalue of the matrix y/p(t)p(t). 

Therefor the analytic expression of the concurrence in 
terms of the parameters defined for the two site system 
is 



C(t) = - 



4ff 



1 + 4.9- 



■) 2 sin 4 (v/TTVJi) . (6) 



We have plotted the evolution of C(t) under RG steps 
versus g for t = 1.5 and t — 0.4 in Fig Q] which shows 
that the concurrence changes from the equilibrium state 
and start to oscillating when the external magnetic filed 
is turned on. As g increases, the height of each peak 
decreases gradually and finally vanish as g — > oo. In- 
creasing the length of chain enhances the oscillation of 
the concurrence versus the magnetic field. However as 
the length of chain increases the first peak of concur- 
rence approaches the critical point (g c — 1) and at the 



3 




01 23456789 10 



FIG. 2: (Color online) Concurrence of the ITF model as a 
function of t for different length chain for g — 0.9. The inset 
shows the concurrence of different lattice sizes collapse on a 
single curve at the critical point (g c — 1). 
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FIG. 3: (Color online) Evolution of the first derivative of T max 
respect to magnetic field (g) under RG. The inset shows the 
Tmax at different RG steps. 



thermodynamic limit the system becomes disentangled 
except at the critical point. Surviving of the concurrence 
at the critical point is the results of the correlation length 
divergence at g c — 1. 

The evolution of concurrence under RG has been plot- 
ted in Fig 12] versus t for g — 0.9 and g c = 1. From 
the Eq.©, it is easy to see that the concurrence is 
periodically fluctuating with the time t with period of 
T = — :22 . Figured shows that the concurrence re- 

duces under RG (increasing the size of system) and dis- 
appears in the large system. But for g c = 1 there is no 
concurrence reduction under RG and the concurrence of 
different length chains coincide with each other (inset of 
Fig. 

The non-analytic behaviour in some physical quantity 
is a feature of second-order quantum phase transition, 
ft is also accompanied by a scaling behaviour since the 
correlation length diverges and there is no characteristic 
length in the system at the critical point. As we have 
stated in the RG approach for ITF model, a large sys- 
tem, i.e. N = 2" +1 , can be effectively describe by two 
sites with the renormalized coupling of in the n-th RG 
step. Thus, the concurrence between the two renormal- 
ized sites represents the entanglement between two parts 
of the system each containing N/2 sites effectively. In 
this respect we can speak of block entanglement -the en- 
tanglement between a block and the rest of system- in a 
large system provided the size of the block and the rest 
of system is equal. 

For any g, there is a time T^ ax (g) at which the C(t) 
reaches its fcth maximum (Fig. (J2J)) and is analyzed as a 
function of coupling g at different RG steps which mani- 
fest the size of system. The first derivative of T^ ax with 

respect to the coupling constant ( dT ™ a * ) shows a singular 



behavior at the critical point. 

We have plotted — for k — 1 (first maximum of 
the C(t)) versus g in FigOfor different RG steps which 
shows the singular behaviour as the size of system be- 
comes large. Surveying the detail shows that the posi- 
tion of the minimum (g m ) of dT ^ g aa: tends towards the 

critical point like g m — gc — N 1 ^ 6 in which 9 = 1 (in- 
set of figure [4]). Moreover, we have derived the scaling 

behavior of In dT ^ g a:l \ gm versus N. This has been plotted 

in Fig.Q, which shows a linear behavior of In \ Qm 

versus ln(N). The scaling behavior is hi dT ™ ax | gm oc N e 
with exponent 9 = 1. 

It is easy to show that the exponent 9 is directly re- 
lated to the correlation length exponent (z/) close to the 
critical point. The correlation length exponent, gives 
the behavior of correlation length in the vicinity of g c , 
i.e., £ ~ (g — g c ) v - Under the RG transformation, Eq. 
([2]), the correlation length scales in the nth RG step as 
C ~ (ffn ~ 9c) ~ v — £,/ n B ' which immediately leads 
to an expression for |^p-| ffc in terms of v and (num- 
ber of sites in each block). Dividing the last equation 
to i ~ {g - g c )- u gives |^| 9o ~ N 1 ^ , which implies 

9 = 1/v, since d " Z ^ ax \ 3m ~ |^rls c a * the critical point. 
It should also be noted that the scaling of the position 
of minimum, g m (inset of figure [3|), also comes from the 
behavior of the correlation length near the critical point. 
As the critical point is approached and in the limit of 
large system size, the correlation length almost covers 
the size of the system, i.e., £ ~ N, and a simple compar- 
ison with £ ~ (g — g c )^ lJ results in the following scaling 
form g m = g c - N 1 '^ . 
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FIG. 4: (Color online) scaling the minimum of — ^ 
various size of system. The inset shows the scaling of the 



of the existences of the finite size scaling for the case of 

block entanglement. We have plotted — d^\g m jg^~ 

versus N{g — g m ) for k = 1 in FigO The upper curves 
which are for large system sizes clearly show that all plots 
fall on each other. 

The similar scaling behaviours and their relation to 
correlation len gth exponent have been reported in our 
previous work [10( in which we have studied the static 
properties of the ground state entanglement of ITF model 
by RG method. 

We would like to mention that in Ref.[l7j the authors 
investigate the dynamics of concurrence of two nearest- 
neighbor sites at ITF model using the exact solution. 

They shows that — f^- dose not diverge at the critical 
point but has a minimum at g = 1. The divergence of 

dT^ 

— at the g c in our work and the similarity of the 
dynamic scaling behaviors to the static scaling behaviors 
f or originates from the low energy state properties. 



position (<?m) of 



dg 



for different length chains. 



IV. SUMMARY AND CONCLUSIONS 



Bp 



"a 




-4 -2 2 4 

N(e-e . J 

■ ' ' ' ' mm 



12 



In this article, we have implement the idea of renor- 
malization group (RG) to study the low energy states 
dynamic of entanglement for spin chains. In this respect 
we show that the RG procedure could be implemented 
to obtain low energy states dynamic of systems in terms 
of effective Hamiltonian which is described by renormal- 
ized coupling constants. This manifest the fact that some 
dynamical quantities of the system could show the finger- 
print of quantum phase transition for an infinite size sys- 
tem. These notions have been observed and approved in 
our study of the ITF model. Moreover, the RG approach 
shows that as the size of the system becomes large, the 
derivative of the time at which the entanglement reaches 
its maximums with respect to the transverse field, di- 
verges at the critical point and its scaling behaviors ver- 
sus the size of the system are as same as the static ground 
state entanglement of the system. 



FIG. 5: (Color online) Finite-size scaling through the RG 
treatment for different lattice sizes. The curves which cor- 
respond to different system sizes clearly collapse on a single 
curve. 



dT 

To obtain the finite-size scaling behavior of — j^Hgmj 
we look for a scaling function in such away that all graphs 
tend to collapse on each other under RG evolution which 
results in a large system. This is also a manifestation 
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